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Abstract 

The newly born millisecond pulsars are investigated as possible energy sources for 
creating ultra-high energy electrons. The transfer of energy from the star rotation 
to high energy electrons takes place through the Landau damping of centrifugally 
driven (via a two stream instability) electrostatic Langmuir waves. Generated in the 
bulk magnetosphere plasma, such waves grow to high amplitudes, and then damp, 
very effectively, on relativistic electrons driving them to even higher energies. We 
show that the rate of transfer of energy is so efficient that no energy losses might 
affect the mechanism of particle acceleration; the electrons might achieve energies 
of the order of 10 18 eV for parameters characteristic of a young star. 
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Introduction 

One of the most fascinating stories in astro-particle physics is the discovery 
of high energy cosmic ray (OR) electrons. Two important constituents of this 
story are: the H.E.S.S (High Energy Stereoscopic System) team’s [1] claim of 
the detection of TeV CR electrons of local origin (within ~ lkpc), and the later 
confirmation of the energy range (.007-lTev) by the Fermi-LAT collaboration 
that analyzed the CR electron/positron data samples collected from 4 August 
2008 to 4 August 2009. [2], 

How do particles acquire such enormous amounts of energy? In astro-particle 
physics, it is conventional to invoke a Fermi-type process: the stochastic ac¬ 
celeration of particles interacting with strong shocks. Though quite effective, 
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such a mechanism requires a pre acceleration of particles to already high en¬ 
ergies [3]. Accelerating light particles like electrons to high energies is even 
more problematic because the strong synchrotron radiation losses will limit 
the maximum attainable energy [4], 

One can imagine electrons being boosted to PeV or higher energies if strong en¬ 
ergy losses could be, somehow, prevented during the accelerating process. One 
such mechanism in which the rotation energy of a neutron star was utilized 
to strongly accelerate particles ( via Langmuir waves generated in magneto¬ 
spheres of pulsars), was recently developed and investigated. It was shown the 
acceleration process proceeds with almost negligible energy losses. [5]. The 
principle steps in the above mentioned process are: 

1) The rotational slow down of a pulsar provides the energy to excite Langmuir 
waves in the bulk electron-positron plasma in the star atmosphere. All pulsars 
are characterized by a decreasing rate of spinning, P = dP/dtt > 0, where 
P is the period of rotation.The period is measured by W = Z0|0|, were 
/ = 2MRI/5 is the moment of inertia of the neutron star, M ~ 1.5 x M 0 
and M 0 ph 2 x 10 33 g are the pulsar’s mass and the solar mass respectively, 

rs 10 6 cm is the pulsar’s radius, = 2n/P is its angular velocity and 

n = dn/dt = -2t tP/p 2 . 

2) The excited Langmuir waves , then, efficiently damp, preferentially on the 
much faster but local beam electrons to accelerate them to even larger energies 

We refer the reader to ([5], and references therein) for the relevant pulsar 
phenomenology, and details of this acceleration mechanism; henceforth, the 
latter twould be termed Langmuir-Landau-Centrifugal Drive (LLCD). LLCD 
is based on the centrifugal acceleration that has been successfully applied to 
some astrophysical settings [6,7]. It is also believed that young pulsars might be 
sources of ultra high energy cosmic rays [8]. For the paper to be self-sufficient, 
we will give a synopsis of the theoretical underpinnings of LLCD before ap¬ 
plying it to the newly born millisecond pulsars; the results are spectacular; 
LLCD can drive electron energies all the way to 10 18 eV. 

In this paper, we do not investigate the aftermath of the acceleration era- 
tor instance, the observational patterns caused by the presence of ultra-high 
energy particles in the magnetospheres of millisecond pulsars. One does expect 
that electrons with such enormous energies will leave their ’’fingerprints” on 
the emission patterns of the corresponding pulsars. The scope of this work 
could be extended to examine the effects of the electromagnetic radiation that 
must come from the ultra relativistic electrons; we intend to deal with these 
problems in near future. Another class of potentially interesting problems ( in 
the context of the present work) are related to the recently discovered events 
of PeV neutrinos by the Ice Cube experiment. Indeed, electrons and positrons 
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with energies much higher than PeV, might produce PeV neutrinos via several 
channels: e~ + e + —> v + e~ + p —» n + v. 

Results 

The Langmuir-Landau-Centrifugal Drive, derived within the framework of a 
relatively simple but nontrivial theoretical model, is shown to work highly ef¬ 
ficiently in the young millisecond pulsars. LLCD, through a two step process, 
converts the pulsar spin-down energy into the kinetic energy of electrons. In 
the first step, the rotation imparted centrifugal particle energy is converted, 
via a parametrically driven ’’two stream instability”, to electrostatic Langmuir 
waves in the bulk electron - positron plasma residing in the pulsar magneto¬ 
sphere. Landau damping of these centrifugally excited electrostatic waves on 
the high energy primary beam electrons (beam) transiting through the same 
region, constitutes the second step. 

The linear growth rates of the rotation driven Langmuir instability, calculated 
for the newly born pulsars, are faster than the typical rotation timescales of 
particles. The Langmuir instability, thus, is efficient and can rapidly convert 
the star slow down energy into the electric held energy 

Landau damping of the unstable Langmuir waves on primary beam electrons 
that converts the electric held energy into particle energy is also shown to 
be equally rapid. The combination creates a very efficient ’’machine” that 
generates ultra high energy (up to EeV) cosmic ray electrons. 

We will show that the LLCD mechanism, when applied to newly born fast 
spinning stars could, indeed, create electron energies that match the highest 
observed in the CR electrons; we believe that this result has much significance 
for high energy astro-particle physics. 

Discussion 

It is generally believed that to attain the possible maximum possible energy, 
the gyroradius of the particle should be contained in the acceleration zone. 
Then, the combination of a strong rotation and a strong magnetic held (of a 
pulsar) leads to an enormous induced electric held on the light cylinder surface 
VtrB/c ~ 10 10 statvolts/cm. This held, might potentially accelerate electrons 
up to energies of the order of E max m 1.3 x 10 19 x Bi c W x P_ 3 eV [ 8 ], where 
Bi c io = Bi c /(10 w Gauss), P _3 = P/(10“ 3 s). What one needs , then, is an 
efficient agency in the rotating magnetosphere that will convert the rotational 
energy into the acceleration of electrons. We believe that the LLCD does 
exactly that: the Landau damping of centrifugally induced Langmuir waves 
guarantees a high efficient channel for acceleration 

For these waves to effectively transfer their energy to particles, the waves phase 
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velocity must be close to the particle speed, which, in this case, approaches the 
speed of light. Further, in the vicinity of v p h, there should be more particles 
a little slower than the wave than particles which are a little faster (in the 
opposite situation, the wave will feed off the particles). For the given problem 
it is always possible to situate v p h in the desired part of the primary beam 
spectrum. Since the distribution function decreases with the Lorentz factor, 
the number of electrons with Vb < v p h exceeds that of the electrons with 
Vb > Vph , where ry denotes the electron speed. Thus the optimum conditions 
for effective Landau damping and, therefore, of net energy transfer from the 
star spin down to electrons, will pertain. 


For relativistic plasmas it has been shown that the Landau damping rate is 
given by [9] 


r ld — 


V'Qj'^b^b 
5 / 2 ~ ’ 

n P l p 


( 1 ) 


where n p is the plasma density and n GJ = BiJfPce ) ~ 1.8 x 10 12 cm 3 is 
the Goldrcich-Julian number density. For the primary beam, 7 b and 07 = 
4ne 2 n GJ /rri are, respectively, the Lorentz factor and the plasma frequency. 
Landau damping is the second step of the acceleration mechanism LLCD. 


We now demonstrate the efficiency of LLCD for young millisecond pulsars 
for which the typical period and slowdown rates are: P ~ 10 _3 s and P ~ 
10 _ 12 ss _1 . It is worth noting that, at this stage such pulsars are only, theoret¬ 
ically, predicted [ 10 ]; this type of pulsars have not yet been observed. 


The efficiency of the LLCD must depend upon the proper coordination of 
the damping rates on the beam electrons and the growth rates of Langmuir 
waves in the bulk plasma. In [5], the magnetospheric plasma was modeled as 
consisting of two streams with the Lorentz factors 7 1j2 (71 < 72 )- It was shown 
that the growth rate of the instability is well approximated by the expression 
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where ay ,2 = y 87re 2 ni i 2/m7 3 2 are the plasma frequencies of the two species, 
n 1,2 are the corresponding number densities, J p (x) is the Bessel function of the 
first kind, b = (2 ck/Q) sin</>_, 2(f)_ = (f> p — (f> e and (f> p and (f> e are, respectively, 
the positron’s and electron’s initial phases. 

It stands to reason that if the instability growth rate was much greater than the 
Landau damping rate, there will be little effective energy transfer from waves 
to the particles. In the diagonally opposite limit with the damping rates far 
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in excess of the growth rates, the waves will not grow much, again resulting 
in very little transfer from the star rotation to the waves. The most optimum 
scenario for an overall efficient energy pumping/transfer system, therefore, 
is realized when the instability growth and Landau damping rates are large 
(with respect to the kinematic rate, hi) and comparable, T ~ T ld- For the two- 
stream instability, the aforementioned condition can be satisfied for particular 
combinations of 71 and 72 ; for example, for 71 « 1.8 x 10 5 , 72 ~ 8 x 10 5 . 

The total energy gained by the beam particles has been estimated to be [5] 

n p F reac 8v 

e ~ — - - 


where Sr ~ c/T, F reac a; 2mcf2£(r) -3 , £(r) = (1 — Q 2 r 2 /c 2 )F 2 [n], From Eq. 
(3) it is clear that two streams with = 7 , fa 1.8 x 10 5 , 72 ~ 8 x 10 5 , which 
in turn guarantee the condition, T ~ T LD , lead to the acceleration of electrons 
up to energies of the order of 10 18 eV. This estimate invokes equipartition 
of energy in the three relevant sp constituents: 7771 fa n 2 72 ~ n GJ ^b with 
7 b fa 7.5 x 10 7 , where n GJ = B/(Pce) fa 1.8 x 10 12 cm -3 is the Goldreich- 
Julian number density [12] in regions in the vicinity of the light cylinder. 

Because the highly relativistic particles will, inevitably, loose energy due to ra¬ 
diation, one must investigate how radiation will affect overall energy transfer. 
Could radiation, for example, put a stringent limit on the maximum energy 
acquired? 

For highly relativistic electrons and photons with ee p h./ (m 2 c 4 ) > 1 (e p h is pho¬ 
ton energy), the inverse Compton mechanism operates in the Klein-Nishina 
(KN) regime [13]. Energy emitted per particle per second is, then, given by 
the approximate expression Pkn oc nrlrn 2 c 5 n p h{e p h) \ In (4ee/m 2 c 4 ) — 11/6 | 
[14] where r e is the classical electron radius. The corresponding time scale 
of the process is Ikn ~ c/Pkn- Using the approximate photon number den¬ 
sity, n p h{e p h ) ~ L/4:7rR 2 c ce p h, and taking into account the typical values of 
luminosity in the high energy level, e p h rs./ lOGeV, satisfy L > 10 35 erg/s, it 
is straightforward to show that for energies of our interest e > IPeV, the 
aforementioned timescale exceeds the instability time scale by many orders of 
magnitude. Thus Compton cooling is very slow compared to the wave energy 
transfer time. The KN time scale goes up with the particle energy implying 
that at higher energies, the inverse Compton mechanism in the KN channel 
is too slow to impose any constraints on the maximum attainable electron 
energy. 

The next possible loss mechanism for relativistic particles, moving in magnetic 
held, is the synchrotron emission with the corresponding estimated power 
P syn ~ 2 e 2 o;^ 7 2 / 3 c [13], where c 03 = eB/(mc) is the cyclotron frequency. 
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According to the standard model of pulsar magnetospheres, there is a signifi¬ 
cant region over the star’s surface, where the electric held is nonzero (vacuum 
gap). It is this electric held that, hrst, accelerates electrons to relativistic 
energies [15]. One can readily show that the electrons leaving the gap with 
a 7 ~ 10 6 will undergo efficient synchrotron cooling at a short timescale, 
t syn ~ 7 me 2 /P syn = 10” 21 s. Thus, immediately after leaving the gap, the par¬ 
ticles radiate their transverse momentum, and very soon transit to the ground 
Landau state. After that, zipping only along the held lines, the electrons will 
reach the light cylinder zone in due course of time. It is precisely the region, 
where the Langmuir waves, always propagating along the held lines, are ex¬ 
cited. Therefore the wave interaction with particles will not cause pitch angle 
scattering, efficiently suppressing the synchrotron mechanism. Quasi linear 
diffusion, another possible source for imparting a pitch angle, [16], also does 
not operate because the required condition for diffusion u> B > uj lj2 , is violated 
for extremely energetic plasmas (for plasmas with energy density exceeding 
that of magnetic held), the synchrotron mechanism is not expected to impose 
any constraints on the maximum attainable energies. 

How about the curvature radiation emitted by particles moving along the 
curved magnetic held lines? To explain force free regime of particles leaving a 
pulsar’s magnetosphere, we have tried to reconstruct the structure of magnetic 
held close to the light cylinder surface [17]. We showed that shown that the 
curvature driven current imparts a toroidal component to the magnetic held. 
As a result, the held lines lag behind the rotation, gradually erasing the insta¬ 
bility. It has been shown that the corresponding timescale is approximately 
given by 


( 2 lb |fcj c 
\3 cu 2 | fcj u b 



k^u b 

~4ff 


Jo 




( 4 ) 


where u b = 7 b c 2 / (co B R B ) is the curvature drift velocity, R B ~ Ri c is the cur¬ 
vature radius of a magnetic held line and and k ± are the components of 
the wave vector of the induced curvature drift mode respectively along the 
drift and perpendicular to it. One can see that for a wide range of param¬ 
eters corresponding to the drift mode with small inclination angles a ( the 
angle of the wave vector with respect to the drift direction), the instability 
timescale is less than the kinematic timescale P. In particular, for a ~ 0.2 
and Ri c k ± ~ 0.1, r CD jrki n <C 1. Therefore, the drift instability is so efficient 
that the magnetic held lines very rapidly reconstruct resulting in a configura¬ 
tion, that enables particles to follow straight line trajectories. The curvature 
radiation is minimized and does not quite interfere with our energy transfer 
mechanism. 

Methods 
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Fig. 1. The distribution function versus the Lorentz factor. It is evident that the 
function consists of two major parts: relatively narrower region with higher Lorentz 
factors, characterizes the primary Goldreich-Julian beam electrons and the wider 
part desc ribes the distribution of electrons resulting from the cascade processes of 
pair creation. 

The LLCD was developed using the standard model of the pulsar magne¬ 
tosphere, where the particle distribution can be, schematically, represented 
as shown in Fig. 1. The narrower region on the figure describes the primary 
beam electrons with higher Lorentz factors ( 7 ?,). The relatively wider sector 
with Lorentz factors in the interval [1; 7 *] (y t < 7 b) characterizes the secondary 
particles produced by means of pair creation. Such a shape of the distribu¬ 
tion function makes the development of two stream instabilities feasible in 
the magnetospheres of pulsars [18]. This feature of the electron distributions, 
along with interesting orbit effects, was exploited in [5] to demonstrate the the 
existence of a strong two-stream like instability that is parametrically driven 
to excite high levels of Langmuir waves . 


In the framework of 1 + 1 formalism [19], the basic system of the Euler equation 
and continuity equations (for each species) and the Poisson equation, is ad¬ 
equate in describing the centrifugally-excited (rotationally driven) Langmuir 
waves. The Fourier transformed linear system consists of 


d P P 

dt 


+ P, 


v P0 tt 2 r B p B + —_E, 




m 


dn p 

dt 


+ ikvpo n p 


+ikn BO v B 


0 


(5) 

( 6 ) 


ikE = 4t vY, n po e p, 

P 


(7) 


where /3 is the species index (electrons and positrons), p B and v p are, respec¬ 
tively, the first order dimensionless momentum (p B —> p B /m) and the zeroth 
order velocity, m is the electron mass, c is the speed of light, is the charge 
of the corresponding particle, ry is the radial coordinate of the corresponding 
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Fig. 2. Here we plot the dependence of k on 72 for three different values of 71 : 
71 = 10 5 (dashed-dotted line); 71 = 1.5 x 10 5 (dashed line) and 71 = 2 x 10 5 (solid 
line). As it is clear from the graph, for a quite wide range of Lorentz factors the 
instability timescale is less than the kinematic timescale, indicating high efficiency 
of the process. 

specie, n 0 and n 0o are respectively the perturbed and nonperturbed Fourier 
components of the density, k is the wave number and E is the electrostatic 
held. The first term in the righthand side of the Euler equation comes from 
the so-called centrifugal force. The centrifugal force, controlled by the follow¬ 
ing ”orbits” of initially relativistic particles r p (t) ~ sin {Vlt + 

v of} (t) ~ V o0 cos (fit + (j)^j [5], where (f) denotes the initial phase of a given 
species, is time dependent. The resulting Langmuir mode equation is peri¬ 
odic in time (has the generalized Mathieu/Hill form) and exhibits parametric 
instability. 

The time dependent centrifugal force that parametrically drives the electro¬ 
static waves is different for the two species - so are their Lorentz factors. 

Since fi r es , the order of the Bessel function is rather high, non-zero growth 
rates will pertain only if the argument b is comparable to fi r es [ 20 ]. 

For the instability to be called ’’efficient”, the inverse growth time Tj ns — 1/Y 
should not exceed the kinematic timescale (also called the escape timescale), 
Tki n ~ P- In Fig. 2, we plot the ratio k = T ins /Tki n versus 72 for different 
values of 7 ^ 7 ! = 10 5 (dashed-dotted line); 7 ! = 1.5 x 10 5 (dashed line) and 
71 = 2 x 10 5 (solid line). The results are derived for typical young millisecond 
pulsars with P ~ 10“ 3 s and Bi c k F* x where B± 10 12 G is the 

magnetic held close to the pulsar’s surface and Ri c = c/VL is the light cylinder 
radius. P 10 _ 12 ss _1 . The condition r ins < r kin is satisfied for a broad range 
of Lorentz factors; the instability is, indeed, efficient. 

Thus we have reviewed, for the young millisecond pulsars, the first step of 
LLCD: the existence of an efficient centrifugally excited parametric instability 
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that pumps the spin down energy of the rotator into electrostatic Langmuir 
waves. Now we review how, through landau damping, this energy will be con¬ 
veyed to the fastest of the particles. The Landau damping in pulsar magneto¬ 
spheres has been investigated in [21], but the authors focus on Alfven waves, 
but not on the electrostatic waves. 
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